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We quantify the dynamical implications of the small-world phenomenon. We consider the generic 
synchronization of oscillator networks of arbitrary topology, and link the linear stability of the 
synchronous state to an algebraic condition of the Laplacian of the graph. We show numerically 
that the addition of random shortcuts produces improved network synchronizability. Further, we 
use a perturbation analysis to place the synchronization threshold in relation to the boundaries of 
the small-world region. Our results also show that small- worlds synchronize as efficiently as random 
graphs and hypercubes, and more so than standard constructive graphs. 

PACS numbers: 



Recently, Watts and Strogatz jjj showed that the ad- 
dition of a few long-range shortcuts to an otherwise lo- 
cally connected lattice (the "pristine world") produces 
a sharp reduction of the average distance between arbi- 
trary nodes. The ensuing semi-random lattice was de- 
noted a small-world (SW) because the sudden appear- 
ance of short paths occurs early on, while the system is 
still relatively localized. This concept has wide appeal: 
the SW property seems to be a quantifiable characteris- 
tic of many real- world structures jl], |[ ||, [Ij , both human 
generated (social networks, WWW, power grid), or of 
biological origin (neural and biochemical networks). 

A spur of ongoing research |2| has concentrated on 
static and combinatoric properties & 3[_ f[ ^, ||, [l^] of a 
tractable SW model |lj [llj . Monasson |11| considered the 



SW effect on the distribution of eigenvalues of the con- 
nectivity matrix (the graph Laplacian) which specifies 
the coupling between nodes — a relevant topic for poly- 
mer networks Jl2|. However, despite their central role in 
real- world networks, there are fewer studies of dynami- 
cal processes taking place on SW lattices. Among those, 
automata epidemics simulations |T^| and Web-browsing 
studies Jl4| have revealed the importance of shortcuts. 
Numerical work on synchronization of Kuramoto oscil- 
lators 0], discrete maps |l5| and Hodgkin-Huxley neu- 
rons |lq| has shown improved SW synchronizability, as 
intuitively expected. However, these numerical examples 
are not generic, and fail to provide insight into how the 
SW property influences the dynamics. 

In this paper, we explicitly link the SW addition of 
random shortcuts to the synchronization of networks of 
coupled dynamical systems. This is an example of dy- 
namics on networks — leaving aside the distinct problem 
of evolution of networks here. By using a generic syn- 
chronization formulation jl7|, [l8| to factor out the con- 
nectivity of the network, we identify the synchroniza- 
tion threshold with an algebraic condition of the graph 
Laplacian. Through numerics and analysis, we quantify 




FIG. 1: Four typical master stability functions (scaled for 
clearer visualization) for Rossler systems: chaotic (bold) and 
periodic (regular lines); with y-coupling (dashed) and x- 
coupling (solid lines) . Here we consider the a;-coupled chaotic 
case (solid bold) with a negative well between (01,02). 



how the SW scheme improves the synchronizability of the 
pristine world, mainly as a result of the steep increase of 
the first-non-zero eigenvalue (FNZE). The synchroniza- 
tion threshold is found to lie in the SW region |j, |l3| , but 
does not coincide with its onset — it can in fact be linked 
to the effective randomization that ends SW. Within this 
framework, we show that the synchronization efficiency 
of semi-random SW networks is higher than standard de- 
terministic graphs, and comparable to both fully random 
and ideal constructive graphs. 

Consider n identical dynamical systems (placed at the 
nodes of a graph) that are linearly and symmetrically 
coupled (as represented by the edges of the undirected 
graph) with global coupling strength a. The topology of 
the graph can be encoded in the Laplacian matrix G, a 
symmetric matrix with zero row-sum and real spectrum 



2 



{Ok}, k = 0, 1, . . . ,n — 1. A general linear stability cri- 
terion for the synchronized state of the system |l7], 
is given by the negativity of the master stability func- 
tion A max (er6>fc) < 0, Vfc. This A max is a characteristic 
of the particular dynamics at the nodes but, crucially, 
a large class of oscillatory systems (chaotic, periodic 
and quasiperiodic) have master stability functions with 
generic features |l8| . In particular, for several chaotic 
systems A max has a single deep well, as depicted in Fig. [l]. 
(We remark that this analysis is quite general: it can be 
extended [ [l8| to eliminate the zero row-sum constraint, 
and to comprise nonlinear coupling and more general syn- 
chronization criteria.) Stability is thus ensured by tun- 
ing the coupling a to try and place the entire spectrum 
of transverse eigenvalues (times a) in the deep, stable 
region: a 9k G (0:1,0:2). This leads to an algebraic con- 
dition for the existence of a linearly stable synchronous 
state: a network is synchronizablc if 



i /9 1 < 0:2/01 = P, 



(1) 



where 9\ is the FNZE and (9 ma x is the maximum eigen- 
value of the Laplacian G. The figure of merit (/?) ranges 
from 5 to 100 for a variety of oscillators (e.g., Lorenz, 
Rossler, double scroll). 

Small-worlds are generated from a pristine world: a 
fc-cycle of n nodes and range k, each node coupled to its 
2 k nearest neighbors for a total of nk edges 11 . The 



Laplacian of this graph G° is a banded circulant matrix 
with non-zero elements on the main diagonal and the 2 k 
adjacent diagonals: G% = 2k and G? = G^ = — 1 with 
(i + 1) mod n < j < (i + k) mod n, and 1 < i < n. 
The SW scheme dopes the pristine world by adding ns 
edges picked at random from the n(n — 2k — l)/2 re- 
maining pairs. Each new edge between nodes I and m 
adds off-diagonal AG; m = AG m / = — 1 and on-diagonal 
AG;; = AG mm = 1 contributions to the Laplacian, thus 
preserving the null row-sum and the bidirectional cou- 
pling. The average number of shortcuts per node (s) is 
related to other measures of randomness (p and q) used 
previously @, 0: s = kp = q(n -2k- l)/(2n). 

The numerical results in Fig. || illustrate the SW ef- 
fect on the synchronization of networks of different size 
and range. For concreteness, all our numerics have been 
performed for a network of identical a;-coupled Rossler 
chaotic oscillators with (3 ~ 37.85. Similarly to other 
locally connected networks, pristine worlds have a large 
eigenratio 6* max /6*i (i.e., they are difficult to synchronize). 
However, as s is increased the eigenratio falls sharply un- 
til, at a value s sync , the condition (|l|) is reached (i.e., the 
addition of shortcuts makes it synchronizablc). The de- 
pendence of s sync on the network parameters {n, k} is 
notably complicated. First, there appears to be an op- 
timal range k ~ 4 for which the SW is most efficient. 
Moreover, the synchronization threshold s sync lies in the 
small- world region but does not seem to coincide with its 
onset. The SW onset (s L ) is defined [Q, by the decay of 
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FIG. 2: Synchronizability thresholds s sync (o) for graphs with 
n nodes in = 300,400,500,1000) and range k £ [1,70], nu- 
merically averaged over 1000 realizations. The solid lines cor- 
respond to the analytical Eq. (^), valid in the range n 1//3 < 
k < fc m i n (n). For most parameters, s sync lies within the small- 
world region between the dashed lines (sl < s < sc) de- 
picted here for n = 1000, but it is distinct from the SW 
onset sl- Note how synchronization is achievable without 
random shortcuts by increasing the deterministic range up to 
fcmm(ro) (see Fig. ^). Inset: decay of the average distance L, 
clustering C, and eigenratio (squares) as shortcuts are added 
to a pristine world of range k = 20 and n = 500. We de- 
fine sl and sc as the points where L and C are 75% of the 
pristine world value; s sync is the point where the eigenratio 
0max/01 =P = 37.85. 



the average graph dista nce L (s) ~ (n/k) /(ns), where 
f(x) = argtanh {xj\Jx 2 + 2x) / \j 'A{x 2 + 2x). Fixing 
L(s L )/L(0) = 3/4, we obtain s L ~ 1.061/n, nk > 1. 
The end of the SW region (sc) corresponds to the effec- 
tive graph randomization through the loss of transitiv- 
ity iQ, ||, ft9||, as given by the decay of the clustering co- 
efficient C|, |o): G(s)/G(0) ~ (2fc-l)/(2fc(l + s//fc) 2 - 
1), n > 1. Again, fixing G(s c )/G(0) = 3/4, we obtain 
s c ^ k (-1 + [{8k - l)/(6fc - 3)]- 1 / 2 ) ~ 0.155 k. The 
synchronization threshold generally lies between these 
two boundaries which scale differently with n and k: 
s L — 1.061/n < s sync < sc — 0.155 k. 

We can gain insight into the synchronization thresh- 
old through an analytical perturbation of the eigenra- 
tio of the SW Laplacian G = G° + G r . Here, G° is 
the deterministic Laplacian of the pristine world, and 
G r is the stochastic Laplacian for the random shortcuts: 
GT. = G r }i = -fy (for i + k + 1 < j < min{n,n-k + i-l} 
with 1 < i < n — k + 1); G^ = (otherwise); and 
Gl = -ELiG« (for 1 < i < n). The are 
n(n — 2k — l)/2 i.i.d. Bernoulli random variables which 
take the value 1 with probability q/n = 2s/(n — 2k — 1) 
(and the value with probability 1 — q/n). The cir- 
culant G° is Fourier-diagonalizable |ll[] with spectrum 
sin[(2/c+l)7rj'/n]/ sin[7r//n], 1 < j < n— 1, 



0? = (2fc+l) 
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(plus 9q = of any Laplacian). The FNZE and the max- 
imum eigenvalue of the unperturbed lattice are: 



9° 

max 



2n 2 k(k 
{2k + 1) -+ 
(2k + 1)[1 



l)(2fc + l)/(3n 2 ) 
3tt/2 



2k- 

•2/37T], 



1 



fc> 1, 



fc<n (2) 
(3) 
(4) 



where ([}]) follows from a continuum approximation. 

Following an "honest" treatment |2lJ with G r as the 
perturbation, we treat the analytical expressions of the 
doubly degenerate eigenvalues as random variables to ob- 
tain their expectations. We postpone the detailed calcu- 
lations |22| and sketch here the main results. After some 
stochastic calculus, the expectations of the eigenvalues of 
the SW Laplacian to second order are shown to be: 



£9 ( P ~ q ± ^3irq/4n 



£9 



(2) 



(5) 
(6) 



for g/n and k/n small. To improve the accuracy of Ssyno 
we have also obtained an approximation to (§) for FNZE: 



S9^ =1 

tUl ~ K 3 



-2q 



~ 9n TS2 f 7 
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(7) 



where if = 2fc+l. Eqns. (0) , (§) , (§ , (|) , and (@) are then 
used to obtain an estimate of s sync as the solution of an 
algebraic equation involving only n and k: 



+f4 1 L = /3(e 1 +M 1) +^i 2) 



(8) 



As shown in Fig. |^, this approximates well our numer- 
ics for n 1 / 3 < k <C n, where the Rayleigh-Schrodinger 
perturbation expansion is valid. 

Using (||), we can obtain 0] a first order estimate 
for the maximum of the synchronization threshold s| ync 
in the valid range. The maximum occurs at k* ~ 
n V (1 + 2/37r)/27r 2 /3 with the asymptotic value s* c ~ 
(2 + 4/3tt)(1 - 2fc*/n)/3(/3 - 1). Therefore, s sync < 
s c [2(l+2/37r)/(2 v / 3-3)( i 9-l)] is bounded by the end of 
the SW region but linked to it. For k < n 1 / 3 , the eigen- 
value bunching (quasi-degeneracy) in the pristine lattice 
renders the doubly degenerate perturbation invalid. We 
are developing another approximation to quantify the be- 
havior in this limit, but our numerics p2j indicate that 
the dependence of s sync with n is sub-logarithmic. This 
confirms that the synchronizability is most effectively im- 
proved for small-range networks (Fig. ||). 

How efficient is the addition of random shortcuts for 
synchronization? We have compared the semi-random 
SW approach with purely random and purely determin- 
istic schemes. An example of the latter is the synchro- 
nization of pristine worlds through the increase of the 
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FIG. 3: Eigenratio decay in an = 100 lattice as / edges are 
added following purely deterministic, semi-random (SW), and 
purely random schemes. Networks become synchronizable be- 
low the dashed horizontal line (/3). The squares (numerical) 
and the solid line (Eq. (|^)) show the decrease of the eigenratio 
of pristine worlds (fc-cycles) through the deterministic addi- 
tion of short-range connections — for n = 100, networks with 
k > 7 are synchronizable. The semi-random SW approach 
(dots, shown for k — 1, 2, 4, 6, 10, 14) is more efficient in pro- 
ducing synchronization. The dot-dashed line corresponds to 
purely random graphs (RG, Eq. ([lo|)), which become almost 
surely disconnected at / ~ 2 In n/(n + 2 Inn) = 0.0843 (thus, 
with 9i G = and unsynchronizable) . The merging of the S W 
and RG behaviors as / — > 1 is the dynamical analogue of the 
effective randomization that leads to sc- 



range k. From (||) and ([|), the eigenratio of a pristine 
lattice is 



9° 

max 



3tt + 2 n 2 



2tt 3 k(k + l) 



(9) 



Therefore, n nodes can be synchronized in a fc-cycle if 
k > fc min ~ ny(37r + 2)/27r 3 /3. (Note in Fig. | the con- 
sistency of our analytical approximation (g): s sync = 
precisely at &; m i n .) For purely random graphs Q n j |p3| , 



9™ x n/ - V2/(l-/)nln7 



9? G nf + V2/F7P" 



(10) 



where / is the number of edges measured as a fraction 
of the complete graph. These are compared with SW 
graphs in Fig. ||. 

We remark on several observations regarding Fig. |^. 
First, the SW addition of shortcuts is more efficient than 
the deterministic addition of short-range layers. Second, 
the effective randomization of SW lattices with edge addi- 
tion translates into converging synchronization behaviors 
of random and SW networks at large /. The / — > 1 region 
is thus robustly stable: cutting connections from the fully 
connected graph has very little effect on synchronization 
stability — not until over 90% are cut (for k small) does 
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connections between topology and dynamics on net- 
works. With a view to improved design, the dy- 
namic eigenratio criterion can be related to other graph- 
theoretical properties (e.g., connectivity, diameter, and 
convergence of Markov chains) [^6). Moreover, other 
measures of cost (e.g., robustness under edge deletion) 
should be considered as possible design constraints. Fi- 
nally, recent results could lead to extensions of this work 
to incorporate more general concepts of stability |27j , and 
broader definitions of small- world lattices jl9| . 

We thank Steve Strogatz for his deep and insightful 
involvement in this work, and Mark Newman for sharing 
computer code and unpublished results. 
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FIG. 4: "Cost" of synchronization measured as the number of 
edges needed to synchronize a lattice of n chaotic Rossler sys- 
tems arranged in different topologies: deterministic graphs (k- 
wheels, fc-cycles, fc-M6bius ladders, bipartites, hypercubes), 
random graphs, and small- worlds (♦). Small- worlds scale fa- 
vorably compared to deterministic structures (and compara- 
bly to the ideal and largely unrealizable hypercubes). Also, 
SW graphs with small range k are as cost-efficient as random 
graphs but demand less (algorithmic) storage memory. 



the eigenratio begin to change drastically. Moreover, if 
we interpret the number of edges needed to synchronize 
n nodes as a simple measure of "cost" , adding many con- 
nections buys little extra stability beyond the small-world 
regime. Finally, it can be shown [£2| that the general 
trends of the eigenratio in Fig. || (namely, "hyperbolic" 
dependence for / small, and near independence for / 
large) can be predicted with the naive perturbation result 
that both 9i and # max change linearly with the number 
of added connections. 

We have also compared the synchronization cost (in 
edges) for SW systems and regular (constructive) lat- 
tices (Fig. |j). For a;- Rossler systems in a fc-cycle, this 
cost tends to / = 2fc m ; n /(n — 1) ~ 0.140. Other con- 
structive lattices |^2|, ^4| also tend to constant fractions: 
/ = 0.252 (for fc-wheels), / = 0.070 (for fc-M6bius lad- 
ders), and / = 0.053 (for the most economical bipartite 
graph). In all those cases, the cost of synchronization 
is high: the necessary number of edges scales like ~ n 2 , 
just like the complete graph. At the other end of deter- 
ministic graphs lies the quasi-optimal (though virtually 
unrealizable) hypercube, which is always synchronizable 
with a number of edges / ~ log 2 n/n. This behavior is 
similar to that of random graphs: from Eq. (|l^) almost 
sure synchronization of Q n j is asymptotically achieved 
when / ~ In n/n. Remarkably, Fig. ^ shows that the cost 
of synchronizing small- /c SW networks is low, i.e., com- 
parable to that of random graphs and hypercubes [p5[ . 

These results hint at research that could deepen the 
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